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Abstract. We resume former discussions of the question, whether 
the spin-spin repulsion and the gravitational attraction of two aligned 
black holes can balance each other. Based on the solution of 
a boundary problem for disconnected (Killing) horizons and the 
resulting violation of characteristic black hole properties, we present 
a non-existence proof for the equilibrium configuration in question. 
From a mathematical point of view, this result is a further example 
for the efficiency of the inverse ("scattering") method in non-linear 
theories. 



1. Introduction 

The examination of time-independent two-body systems dates back to the early 
days of General Relativity. In a 1922 paper, Rudolf Bach and Hermann 
Weyl J3] discussed the superposition of two exterior Schwarzschild solutions in 
Weyl coordinates as a characteristic example for an equilibrium configuration 
consisting of two "sphere-like" bodies at rest. Bach noted that this static solution 
develops a singularity on the portion of the symmetry axis between the two 
bodies, which violates the elementary flatness on this interval. In a supplement to 
Bach's contribution, Weyl focused on the interpretation of this type of singularity 
and used stress components of the energy-momentum tensor to define a non- 
gravitational repulsion between the bodies which compensates the gravitational 
attraction. Weyl's result is based on some artificial assumptions but implies an 




Figure 1. Illustration of two aligned rotating black holes with horizon areas A\, 
A2 and angular momenta J\, J2. 

interesting question: Are there repulsive effects of gravitational origin which 
could counterbalance the omnipresent mass attraction? 

Post-Newtonian approximations tell us that the interaction of the angular 
momenta of rotating bodies ("spin-spin interaction") could indeed generate 
repulsive effects. This is a good motivation to study, in a rigorous way, stationary 
two-body problems. 

In this contribution we shall summarize the results that we obtained for a 
stationary two-black-hole system consisting of two aligned rotating black holes 
with parallel (or anti-parallel) spins, see Fig. Q] The representation is based 
on three recent papers which contain the details of the analysis ||29] [T4l [30l . 
The idea of our non-existence proof is to construct the exterior gravitational 
fields of two disconnected Killing horizons, see Fig. [2] below, via a boundary 
problem for the nonlinear Ernst equation, which is essentially equivalent to 
the vacuum Einstein equations. Fortunately, this equation belongs to a class 
of completely integrable differential equations, which can be mapped to linear 
structures ("Linear Problems"). This fact is the source of powerful solution 
generating methods such as Backlund transformations. It can be shown that 
a single Backlund transformation ||251 1271 applied to Minkowski space creates 
a Kerr-NUT spacetime which includes the spacetime of the rotating black 
hole. Since iterative Backlund transformations act as a "nonlinear superposition 
principle", the double-Kerr-NUT solution [19. 25 1 was considered to be a good 



candidate for the solution of the two-horizon problem and extensively discussed 
in the literature O] [J5] [16] H7] US] [J9] |20] EI] |22] [37] SO]. However, there 
was no argument ensuring that this particular solution be the only candidate. 
We have removed this objection and shown that the solution procedure for 
the boundary problem necessarily leads to a subclass of the double-Kerr-NUT 
solution. Thus we could make use of results derived for the double-Kerr-NUT 
solution. The result is in line with a theorem of Varzugin ll38l [39], which 
says that the 2A^-soliton solution by Belinskii and Zakharov ||U [5) contains all 
possible solutions (if any existed) corresponding to an equilibrium configuration 
of black holes. The subclass is characterized by a set of restrictions for the 
parameters of the general double-Kerr-NUT solution. These restrictions, first 
derived and discussed by Tomimatsu and Kihara 111711371 . ensure the regularity of 
the double-Kerr-NUT solution on the axis of symmetry and on the horizons. An 
elegant reformulation of the Tomimatsu-Kihara regularity conditions by Manko, 
Ruiz and Sanabria-Gomez [22] made it possible to express black hole quantities 
such as mass, angular momentum and surface area in terms of independent 
parameters. We have made use of these results. Another condition to be satisfied 
is the positivity of the total mass. Combining the restrictions with symmetry 
arguments, Hoenselaers and Dietz ifTT] [T6l and Krenzer 11201 could show that 
the double-Kerr-NUT solution cannot describe a configuration consisting of two 
identical black holes. Manko and Ruiz 1211 generalized this result by showing 
that the regularity conditions imply that at least one of the two horizons has a 
negative Komar mass. They argued, without giving any explanation, that this 
peculiarity casts out the double-Kerr-NUT solution. Remarkably, Ansorg and 
Petroff [1], who described an equilibrium configuration with a positive total 
mass and a component which has a negative Komar mass, came to an opposite 
interpretation. However, considerations like these stimulated us to examine 
further black hole inequalities. Fundamental local "state variables" of a rotating 
black hole are its area A and its angular momentum J. Indeed, for a single 
black hole these quantities are restricted by the inequality 8n\J\ < A. Based on 
results of Ansorg and Pfister [2], who examined extremal black holes, Hennig, 
Cederbaum and Ansorg [13], who, following Booth and Fairhurst [6|, studied 
sub-extremal black holes defined through existence of trapped surfaces (surfaces 
with a negative expansion of outgoing null rays) in every sufficiently small 
interior neighbourhood of the event horizon, and Chrusciel, Eckstein, Nguyen 
and Szybka [8 |, we can assume that each of the two black holes has to satisfy 
the angular momentum-area inequality individually. Surprisingly, Dain and 
Reiris [10] were able to extend its range of application to non-stationary black 
holes, see also the overview article by Dain [9] and references therein. 



2. Mathematical tools 

2.1. Metric and horizons 

The exterior vacuum gravitational field of axially symmetric and stationary 
gravitational sources can be described in cylindrical Weyl-Lewis-Papapetrou 
coordinates (gX,<p,t in which the line element takes the form 

ds 2 = e- zu [e 2k (dg 2 + dC z ) + Q 2 dcp z ] -e 2u {dt + adtp) 2 , (1) 

where the "Newtonian" gravitational potential U, the gravitomagnetic potential 
a and t he "supe rpotential" k are functions of g and £ alone. At large distances 
r = I V 6 2 + C 2 1 —I °° from isolated sources located around the origin of the 
coordinate system, r = 0, the spacetime has to be Minkowskian, 

r -)• oo : di 2 = dg 2 + dC 2 + g 2 d<^ 2 - df 2 . (2) 

Metric (fTJ admits an Abelian group of motions G2 with the generators (Killing 
vectors) 

C = 6' t , (stationarity), (3) 
rf =8', (axisymmetry), (4) 

where the Kronecker symbols 5), 5 l indicate that has only a time f-component 
whereas rf points in the azimuthal ((/?-) direction, rf has closed compact 
trajectories about the axis of symmetry and is therefore spacelike off the axis 
(and the horizons). £' is timelike sufficiently far from the black holes but can 
become spacelike inside ergoregions. Obviously, 

e 2u = a = - e - 2U Vi C (5) 

is a coordinate-free representation of the two relativistic gravitational potentials 
U and a. 

In stationary (and axisymmetric) spacetimes, the event horizon as the central 
black hole property is a local concept. Consider the Killing vector 

Z' = £ + nri (6) 

with the norm 

where f2 is a real constant. A connected component of the set of points with 
e 2v = 0, which is a null hypersurface, (de 2v ,de 2y ) =0, is called a Killing horizon 

M?(g): e 2y = -(e',e')=0, (de 2y ,de 2v )=0. (8) 

1 In the following, we also use the complex coordinates z = Q + iC and f = ^ — if. f is the 
time coordinate. 
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Figure 2. A two-black-hole equilibrium configuration in Weyl-Lewis- 
Papapetrou coordinates. (Adapted from 11291 .) 



Since the Lie derivative jSf^/ of s 2v vanishes, we have (£',de 2y ) = 0. Being null 
vectors on Jf?(£'), £' and de 2v are proportional to each other, 

(O : de 2v = -2k£. (9) 

Using the (vacuum) field equations one can show that the surface gravity k 
is a constant on Q is the angular velocity of the horizon. In black hole 

thermodynamics, k and Q are conjugate to the extensive quantities A (area) and 
J (angular momentum), respectively. 

In the q-C, plane (f = constant, \p = constant) of the Weyl-Lewis-Papapetrou 
coordinate system (03, the horizons cover a finite portion on the £-axis (q = 0), 
see Fig. |2] or shrink to a point [7 |. It turns out that extended horizons ("sub- 
extremal horizons") and point-like horizons ("degenerate horizons") require 
different considerations. Note that a Killing horizon is always a two-surface in 
the time slice t — constant. The degeneracy to a line or a point is a peculiarity of 
the special coordinate system. 

In this paper, we explain the non-existence proof for extended (sub-extremal) 
horizons and end up with a brief comment on degenerate horizons. 

The dashed line in Fig. [2] sketches the boundaries of the vacuum region: jz/ + , 
are the vacuum parts of the £-axis (axis of symmetry), Jff^ and Jf 1 ^ 1 
denote the two Killing horizons, which are located in the intervals [X^i^i] and 
[^4,^3] on the ("-axis, and ^ stands for spatial infinity. The gravitational fields 



a, k, U have to satisfy the following boundary conditions 



si*,!* u : a = 0, k = Q, (10) 

l + A-a = 0, ? = 1,2, (11) 
<£: U^O, a^O, fe-)-0, (12) 

where fl\ and i?2 are the angular velocities of the two horizons. Equations (1 1 Ob 
characterize the axis of symmetry (rotation axis). The first relation originates 
from the second equation in (|5}, since the compact trajectories of r] with the 
standard periodicity 2tt become infinitesimal circles with the consequence 77 — ^ 0. 
The second relation is a necessary condition for elementary flatness (Lorentzian 
geometry in the vicinity of the rotation axis). Equation (fTTT i is a reformulation of 
Eqs. © (e 2v = 0) and © since the horizons are located on the £-axis (g = 0), 
see Fig. 12 Finally, Eq. ( 112b ensures the asymptotic flatness of the metric (fTJ, see 



2.2. Field equations and Linear Problem 



The stationary and axisymmetric vacuum Einstein equations for the metric 
potentials U and a are equivalent to the Ernst equation 



W)(f.. ee +f.cc + - g f. e ) =f 7 e +f: 

for the complex function 

f(Q,C)=e 2UieX) +ib(e,C), 



where the twist potential b is defined by 



. e = ge 4u b x , 



-oe- 4u b, e . 



The potential k can be calculated from 



(13) 

(14) 
(15) 

(16) 
(17) 

As a consequence of the Ernst equation ( 113b . the integrability conditions = 
a^ e and k, B ^ — k ^ e are satisfied such that the metric potentials a and k may 
be calculated via line integration from the Ernst potential /. Since e 2u = 
all metric coefficients in (HJ can uniquely be determined from /. Thus the 
boundary problem for the vacuum Einstein equations reduces to a boundary 
problem for the Ernst equation. However, we have to cope with non-local 
boundary conditions (fT0li-(fT2l. (fT5ll-(fT7li for the Ernst potential /. Fortunately, 
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these boundary conditions are well-adapted to the "inverse method", which we 
applied to solve the boundary value problem. 

The Ernst equation is the integrability condition # ,j = ^ j z of the Linear 
Problem (LP) 



(18) 
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(19) 



where the pseudopotential £»(z,z,A) is a 2 x 2 matrix depending on the spectral 
parameter 



K-iz 
K + iz 



as well as on the complex coordinates 

Z = Q+i(, 



(20) 



(21) 



whereas M, N and the complex conjugate quantities M, are functions of z, z 
(or g, alone and do not depend on the constant complex parameter K. Since 
the integrability condition must hold identical in A (or K) it yields a system of 
first order differential equations for and M which is equivalent to the Ernst 
equation. The first order system has the "first integrals" 



M = 



f,z 
f + f' 



N = 



h 
f + f 



(22) 



Vice versa, any solution / of the Ernst equation admits the unique determination 
of the pseudopotential # up to constants of integration. Thus the Ernst equation 
( fT3l and the LP ( fT8l ). ( fT9b are equivalent to each other. 

Multiplying ( TT~8T > by dz and adding ( fT9b multiplied by dz one obtains the 
reformulation d$ = (•••)$ of the LP in the form of a system of (overdetermined) 
total differential equations. 

Without loss of generality we choose the standard representation 



x(^,C,A) -x(£>,C.-A) 



where 



■0 e,Cj- 



i 



: x(e,C 5 A) 



(23) 



(24) 



due to the special structure of the coefficient matrices of the LP. For K — > °° and 
A — >• — 1, the functions ip, \ can be normalized by 



^(ftC,-i) = x(e,C,-i) = i. 



(25) 



As a consequence of the LP, the Ernst potential and the gravitomagnetic potential 
can be read off from the pseudopotential # at K — > °° and X — > +1, 

/(e,0=xG?,C,i) = *2iG?.C,i), (26) 

X = 1 ~ c < < 27 > 

where C is an arbitrary constant. The idea of the inverse (scattering) method 
is to discuss ^, for fixed but arbitrary values of g, £ (z, z) as a holomorphic 
function of X or ZT. In the latter case, # "lives" on the two sheets of the Riemann 
surface associated with ( 1201 i. As this mapping depends on the parameters g, £, 
the position of the branch points Kb = if, A"b = — iz and the branch cut between 
them changes with the coordinates. 

3. Non-existence proof 

3.1. Integration of the Linear Problem 

In order to solve the boundary problem (TT3K (ITUi - (TTZt , we shall integrate 
("solve") the LP along the dashed line in Fig. [2] which marks, in the plane, 
the boundary of the vacuum domain outside the horizons. Starting from and 
returning to any axis point, say g = 0, £ 6 &/ + , we shall make use of the boundary 
conditions and finally arrive at a representation of the Ernst potential on the 
axis of symmetry. It turns out that this representation is sufficient to express 
all black hole quantities (such as areas A i, A 2 and angular momenta J\, J2 of the 
black holes) in terms of three independent real parameters (plus two additional 
scaling parameters) and to establish the equations of state of the black hole 
thermodynamics of the equilibrium configuration under discussion. Furthermore, 
the axis values /(£) = f(g = 0,0 fix the solution f(g,() of the Ernst equation 
uniquely |[T2| . 

Since X(K) as defined in (l20b degenerates on the £-axis, X = ±1, the LP 
d^ = (...)^ can easily be integrated. For X = +1 one obtains 

^,.<^>: •_(/ _\)L, L-($> (28, 

The representation for X = — 1 follows from ( 123b by exchanging column 
elements. Remarkably, the ^-matrix separates. The first factor depends on 
the path of integration whereas L representing the "integration constants" is a 
function of the spectral parameter K alone. There is a difference between the 
case of two extended horizons and that of one or two degenerate horizons. In the 
first case one can parametrize the dashed curve by the coordinate £ everywhere on 
the £-axis, i.e. f=f(g = 0X)=f(() on .e^,^ ,^''). This is clearly impossible 
if the dashed curve runs around a point-like horizon. A path like this can be 



-2V 



K 1 



dK 



[ x (-a)-v(-a); 



described by an infinitesimal semicircle which brings (local) polar coordinates 
(R, 9) into play (24]. Then we have f = f(R-> 0, 9) = f{9), 9 6 [0,vr] in ((28). 

To exploit characteristic properties of the horizons such as (© and ©, it is 
helpful to introduce corotating frames of reference defined by 

g' = g, C' = C, tff = <p-nt, t'=t, (29) 

where Q = Q\, J?2- This coordinate transformation induces transformations of 
the gravitational potentials in ((TJ such that 

e 2U> = e 2V = e 2i/ [(1 + Qa) 2 _ Q 2 Q 2 t - W] ^ (3Q) 

{\-nd)er u ' = {\ + na)e 2u , (31) 

where a prime denotes "corotating" quantities. To determine the corotating Ernst 
potential /' one has to apply d 1 5b to a' and e 2u . Finally, using the equations (122b 
for N',M', one obtains the corotating pseudopotential [24 1 



= T n #, (32) 



where 



+ Qa- f2ge- 2u 
Ln ~ { 1 + fla + n g e- 2u 

+i(K + iz)f2e- 2u (^ ^ ~ A V (33) 

The validity of the Ernst equations of the non-rotating and corotating system 
at the points of intersection g/ /Jtf (axis of symmetry/extended horizon) and 
srf 1^) (axis/circle at infinity), see Fig. [2] implies that ^ and must be 
continuous there as well. By way of example let us consider the continuity of 
^ in d28T i and in ( 1321 ) at the point g = 0, £ = K\ . It immediately leads to a 
connection of the horizon and axis values of the "integration constants" LA 1 ) and 
L+, cf. (|28]l, 

L«=fl -W, (34) 

V 2i/?( l )(^- J p 1 )y ' 

where Fi is is a special case of F, as needed later, 

Fi := ( ~j| * V /, = / = 1 , • • • ,4. (35) 

Note that fx = —fx since e 2y and e 2t/ are continuous at the points of intersection 
and e 2v = from the side of the horizon. If one continues the interlinking 
procedure along the closed contour, one returns to the starting point with the 
result 



1 

1 



(L+)- 1 =R+, (36) 



where 

R+ -n( 1 -<-')'5»^)(: i) 

with nO) = i2(2) = n u «(3) = fl(4) = 

Point-like (degenerate) horizons can be involved without any difficulty by 
setting K{ = K% or/and Kj, = K4 in these equations \24\. 

We shall show that ( 136b with ( f37b as the result of the integration of the LP 
along the closed (dashed) contour in Fig.|2]yields the Ernst potential on the axis. 



3.2. Ernst potential on the axis 

At the branch points K = Kb = v., K = Kb = — h. of the Riemann /f-surface, 
where X = 0, A = 00, respectively, one finds from (1231 <&u = &12> &2\ = —&22, 

'-*=•(! i>--(i .,)■ 

On the £-axis one has confluent branch points, Kb = Kb — (. For this choice one 
obtains from ( 1381 l. ( 128b and ( I361 l in terms of the Ernst potential / + on s</ + 

and so 

[R + (C)-i]( f+ \ c) ) =0 (40) 

for R+(C) := R + (K = C), cf. (El). Note that gD) is equivalent to OS: a second 
column resulting from (l39l l and (140) are complex conjugate. 

We shall now discuss properties of the Ernst potential on the axis which can be 
derived from the eigenvalue equation (140) . First of all, let us point out that similar 
equations can be derived for all intervals stf ', Jti? by the interlinking procedure 
as explained in ( 134b . At the first glance, / + (C) seems to be a quotient of two 
polynomials of fourth degree in £. However, attention must be paid to the fact 
that the characteristic determinant has to vanish, 

|R+(C)-1|=0. (41) 

This condition tells us that the numerator and the denominator of / + (C) must 
have two common zeros such that the axis potential is a quotient of two 
(normalized) polynomials of second degree in £, 



where q, r, s, t are complex constants which can be expressed in terms of /,, 
Kj and J?M, (i = 1, . . . ,4). For extended (sub-extremal) horizons the following 
reparametrization is useful: 
Defining 



a: := 



d 2 {Kj) 
d 2 (KiY 



OCiCCi = l, fii\= 



n 2 (Ki, 



and using 



one obtains 



-f + (K,), 
Pi = -oti. 



n 2 {KiY 
i=l,...,4 



m = i 



(43) 

(44) 
(45) 



The equations (l43l form a linear algebraic system for the parameters q, r, s, t. 
Eliminating them in ( 142t one arrives at a determinant representation for the Ernst 
potential on the axis , / + (C), which can be written in the form 
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ij=l,...,4. 



(46) 



(47) 



Rij := 



./'./ 



where 



(C-#,) 2 + ? 2 >0, '../ 1 4. 



(48) 



(49) 



f(g,C) I s tne Ernst potential of the double-Kerr-NUT solution which was 
originally generated by a two-fold Backlund transformation of Minkowski space 
Il25lll9l in the form of a quotient of two 5x5 determinants. According to a rule 
of Yamazaki [40], this type of determinants can be reduced to 2 x 2 determinants 
as used in tgS). Making use of 02) [or (|27)], 02), 02) and e 2u = 91/, one finds 
determinant representations for all metric coefficients, i.e. for a, k, e 2u in 0~K 
see 1301 . 

As a condition identical in K, Eq. (1411 i yields four constraints among the 
parameters Q\, fl 2 ; K\ — K 2 , K 2 ~K^, —K4; f\, . . ., fn,. It can be shown QUI 



that this system of algebraic equations guarantees that a = on si ± , £/° and % '. 
As a consequence, it eliminates NUT parameters from the Ernst potential. In 
consideration of (l40l > and introducing dimensionless coordinates g, £ via 



e=r-> C = hf — (50) 
A23 A23 

one realizes that f + (() from ( 140b and therefore f(gX) depend on four real 
parameters. 



3.3. Weyl-Bach force between the black holes 

So far we have examined the Ernst equation as a classical field equation. It is 
questionable whether the parameter conditions (HTt alone could rule out the Ernst 
potential under discussion and lead to a non-existence proof. Consider the static 
Ernst equation which is an axisymmetric Laplace equation for the "Newtonian" 
gravitational potential U in Weyl coordinates. The superposition of two solutions 
with aligned rod-shaped sources (as "classical" precursors of horizons) is regular 
outside the sources. It is the gravitational interaction ("force") between the 
rods that forbids equilibrium. Bach, who examined this example in the already 
mentioned Bach/Weyl paper J3|, noted that the metric function k cannot vanish 
on the portion of the axis between the two sources and that this fact violates the 
regularity of the solution. Weyl's remarks (published as a supplement to Bach's 
paper) focused on the interpretation of this type of singularity. He used fictitious 
stresses described by an energy-momentum tensor to define a force of attraction 
between the sources. This "Weyl-Bach force" turned out to be proportional to 
a constant value of k on the portion of the axis between the two sources (in our 
notation k°) provided that k + = 0, k~ = 0, which is a possible gauge. Note that k 
as defined in ( II6I 1. (fTTT i has constant values on all intervals of the £-axis. One of 
them, considered to be the arbitrary integration constant, can be chosen so that, 
say, k + = 0. Integration of (fTol i. (fTTT i along ^ in|2]results in k~ = 0. 

Equipped with that physical as well as geometrical interpretation of the 
superpotential k° ("attractive force" that violates the Lorentzian geometry in the 
vicinity of the rotation axis) we shall examine the boundary condition 

k°=Q (51) 

for the gauge 

k+ =k~ = 0. (52) 

Our discussion is based on the original parametrization of the double Kerr- 
NUT solution, see d43}, (|45j, (|46j, and Kramer's representation of e 2k £18]. (In 
principle, one could determine the axis values of this gravitational potential from 
the axis values of the Ernst potential by integrating the equations (fTol i. (fTTT i along 
the axis, i.e. by operations on the dashed contour in Fig. [2]) It turns out that 



the boundary conditions for non-overlapping extended horizons (fSTb . ( 1521 > are 
satisfied by only one parameter condition: 



oi«2 + 0130:4 = 0. (53) 

Two of the four conditions (HTt can be used to eliminate Q\ and J?2- The two 
remaining equations turn out to be equivalent to the equations 



(I-04) 2 2 (I-03) 2 /H^Mer, N ,™ 

w = , w:= 4/ — — g [1,°°), (54) 



04 03 V ^13^23 

(l+o 2 ) 2 /2 _ (l+O!) 2 , /^K^ 



w' 2 = ^ ^, w ':=W-^-^e(0,l]. (55) 

«2 Oi V ^13^14 

The three restrictions (l53~t . (I54K (l55l are nothing else but a reformulation of 
the original Tomimatsu-Kihara conditions IfTTl [371 . This reformulation is due 
to essential examinations of the double-Kerr-NUT solution by Manko, Ruiz and 
Sanabria-Gomez B221 |2D . As was particularly shown in iBTI . the restrictions 
( T53I ). ( I541 i. (IB3T l can be solved to express the parameters a\, . . . , 04 in terms of the 
three real parameters w, w' and 0: 

w'a 2 + iea a 2 + iw'ea 

01 = — ; — : , o 2 = - — — , (56) 

w — lea 1 — w'ea 

wa 2 — a a 2 — wa 

03 = , o 4 = - , (57) 

w — a 1 — wa 

where e = ± 1 and 

0304 = -0102 = o 2 (o = e i0 , (j)£ [0,2tt]). (58) 

Now we have arrived at the final form of the solution of the boundary problem 
( TTOb . ( fTTI ). ( 112b for the Ernst equation (113) . Eliminating the a,- in favour of w, w', 
4> and introducing dimensionless coordinates (l50l l. / becomes a function of two 
coordinates and three real parameters^, 

f = f(~gX;w,w'^). (59) 

Note that the relative horizon "lengths" K12/K23, ^34/^23 can be expressed in 
terms of w,w' as well. At this point we cannot guarantee that this Ernst potential 
is well-behaved. Computer experiments show that the regularity of the Ernst 
potential on the axis of symmetry must be "paid of" in the form of singular rings 
outside the horizons. Fig.[3]conveys an impression of the structure of this type of 



2 From this point of view, the Ernst potential of the Kerr solution depends on one real 
parameter and two dimensionless coordinates. 



Figure 3. The singular Ernst potential for a particular example configuration. 
Parameters: <j> = -0.1, w = 1.3, W = 0.5, K\ = 2, K 2 3 = 2, e = 1. 



singularity. Irregularities on the horizons could also rule out the solution. Since 
area A and angular momentum J are individual quantities of each black hole, it 
is obvious to examine the inequality 87r|/| < A which characterizes a single sub- 
extremal black hole. Based on the literature as commented on in the introduction, 
we can take for granted that this inequality has to hold on both horizons of a 
regular spacetime with two sub-extremal black holes. 

3.4. Angular momentum-area inequality and non-existence proof 
In order to examine the inequalities 

JF®: Sir\Ji\ <A U z = l,2, (60) 
we calculate the quantities 

Pi-=— — , ' = 1,2, (61) 

in terms of the parameters w, w' and (f> € [0,2ir). Note that this computation can 
be performed with the aid of the axis values of the Ernst potential alone, see 1291 . 
The result is 

l + <Pw' w(w-<P) 

px=£ w i{<p +w >y P2 = £ T^r (62) 



where 

From this we have 
and 



<P := cos</> + £siri(/>, e = ±l. (63) 



P f-l = (l-w^) ^^±l <0 (64) 
rl v ; w' 2 (<£ + w') 2 

p|-l = y-l)_^-<0. (65) 

For the allowed parameter ranges w 6 [1,°°), w' G (0,1], which follow from the 
definitions 654\ . (IB3T > of w and w', these inequalities can only hold if 



w 



' 2 1 2<£w' + 1 < and w 2 -2&w + 1 < 0. (66) 



This, however, implies <Pw' < and <Pw > in contradiction to w' > and 
w > 0. Thus we have proved the non-existence of stationary and axisymmetric 
configurations consisting of two aligned sub-extremal black holes and we 
conclude that the spin-spin repulsion cannot compensate for the gravitational 
attraction. 



4. Summary 

As a characteristic example for the ongoing discussion about existence or non- 
existence of stationary equilibrium configurations within the theory of General 
Relativity, we have studied the question whether two aligned, sub-extremal black 
holes can be in equilibrium. The result of our above analysis, whose details can 
be found in 11291 1141 1301 . is negative: there are no two-black-hole equilibrium 
configurations! The idea of the non-existence proof is illustrated in Fig. [4] and 
can be summarized as follows. 

Equilibrium configurations with two aligned rotating black holes, if any 
existed, can be described by a boundary value problem for two separate (Killing-) 
horizons. Remarkably, this problem can be solved by integrating the Linear 
Problem d^ = (...)# along the dashed contour as sketched in Fig. [2] or Fig. [4] 
Thus we arrive necessarily at particular Kerr-NUT solutions which have two 
horizons [at least according to the definitions (0 and ©] and show the correct 
regular behaviour at infinity and on the symmetry axis, whereas regularity off 
the axis is not guaranteed. On the contrary, we find that all candidate solutions 
indeed do suffer from irregularities. One of them is the violation of the angular 
momentum-area inequality 8ir\J\ < A, which must hold for any regular sub- 
extremal black hole. However, we could show that there is no choice of 
parameters for which angular momentum and area of the two horizons jointly 
satisfy the inequality. Hence, there exists no regular solution of the vacuum 
field equations for stationary two-black-hole configurations. For brief comments 



d#= (...)# : + Sn\J\<A 



Non-existence of 
two-black-hole equilibrium 



Figure 4. Summary of the non-existence proof for two sub-extremal black holes. 



on the extension of the non-existence proof to degenerate black holes, see the 
following supplement and Fig. [5^. 

5. Supplement 

5.7. Degenerate horizons 

The analysis as presented above only applies to configurations with two sub- 
extremal ("extended") horizons. But we have already indicated several times 
that it is possible to extend the proof to configurations containing one or even 
two degenerate black holes (with "point-like" horizons, where K\ = K2 or/and 
K3 = K4), see Fig. [5^. In the first part of this supplement we give an outline of 
this generalisation. For details we refer to lTl4l . 

As already mentioned, the degeneracy of a Killing horizon to a point is merely 
a peculiarity of the Weyl-Lewis-Papapetrou coordinates. In order to resolve the 
internal structure of the horizon, we have discussed the LP in polar coordinates 
centred at the point-like horizon in Ifl4l . Integrating the LP, we found that all 
possible equilibrium configurations with degenerate black holes can be obtained 
as particular limits of the double-Kerr-NUT family of solutions. It turned out that 
we obtain two families of (two-parametric) solutions describing configurations 
with one degenerate and one extended horizon, whereas we found three (one- 
parametric) solution families for configurations with two degenerate horizons. 
In order to exclude these families as acceptable equilibrium configurations, we 
showed that they suffer from unphysical singularities. 

In the case of one degenerate and one sub-extremal black hole, the angular 
momentum-area inequality becomes an equality for the degenerate horizon 
(which turned out to be satisfied identically and, therefore, did not provide 
any new information). But it still has to hold as an inequality for the sub- 
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Einstein-Maxwell 
equations 

3x3 Linear Problem 



Figure 5. Illustration of further BVPs that can be solved with the inverse 
scattering method: (a) two-horizon problem with one degenerate black hole, 
(b) rigidly rotating disk of dust, (c) constructive uniqueness proof for the Kerr- 
Newman black hole. 



extremal black hole. This restricts the parameters, but does not yet exclude the 
possibility of a regular equilibrium configuration. Hence an additional ingredient 
was required for the desired non-existence proof, which was the positivity of 
the total mass (ADM mass) of the spacetime. As first shown by Schoen and 
Yau 1351 l36l . the mass of a regular, asymptotically flat spacetime satisfying the 
dominant energy condition (which is certainly satisfied for a black hole vacuum 
spacetime) is positive. However, for the configuration in question, we found 
that the entire parameter range in which the sub-extremal black hole satisfies the 
angular momentum-area inequality has a negative mass. Thus we could conclude 
that the solutions with one degenerate and one sub-extremal horizon are singular. 

Finally, we studied the three solution branches for configurations with two 
degenerate horizons. One branch could be excluded since it has a negative ADM 
mass for all possible parameter values. The other two branches have negative 
masses for most parameter values, but there are small parameter regions with 
positive mass. Fortunately, these solution branches are relatively simple and, by 
studying the solutions of a certain quartic equation, it was possible to demonstrate 
explicitly that singularities (singular rings around the £-axis) are present even in 
the parameter range with positive ADM mass. 

Hence we could extend the non-existence proof to all forms of horizons. 

5.2. Further applications of the inverse method 

In the following we briefly comment on some other applications of the inverse 
(scattering) method to rotating objects. The integration of the Linear Problem of 
the Ernst equation d# = . . . along a suitable closed contour was first practised to 



determine the gravitational field of a rigidly disk of dust 031113211331 . see Fig.|5]3. 
Among other things, the solution (expressed in terms of theta functions) describes 
a parametric collapse of the disk with a final phase transition to an extreme black 
hole. 

Obviously, the integration method under discussion can be used to construct 
the Kerr solution as the unique solution of the one-horizon boundary problem 
ll28l . This corresponds to the methods of electrostatics and is an alternative to 
"complex tricks" and other formal derivations of the Kerr metric in the textbooks. 

Extending the vacuum examinations to Einstein-Maxwell fields, Meinel 
11231 has recently constructed the Kerr-Newman solution by integrating a 3 x 3 
electrovacuum LP, see Fig.[5j: and thus proved the uniqueness of the solution. 
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